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AB ST RAC T 



The number of zeros contained in a given interval (a,b) for a 
solution to a Sturm-Liouville differential equation is of importance in 
many problems of mathematical physics. This number maybe deter- 
mined through Sturm’s Comparison Theorem. Given one zero of the 
solution to a Sturm-Liouville differential equation, a technique, based 
upon Sturm’s Theorem, of computing the next consecutive zero of the 
solution is proposed. The existence of a function which satisfies the 
desired end results of the proposed technique is shown. The technique 
is then applied to Bessel’s differential equation and the results tabula- 
ted for the first 20 roots of and Jj(x), Unfortunately this 

technique did not achieve the desired result of convergence to succes- 
sive zeros of the given Bessel Function, 

The writer wishes to express his appreciation for the assistance 
and encouragement given him by Professor E. J. Stewart of the U. S, 
Naval Postgraduate School in this investigation. 
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1. Introduction, 






Consider a differential equation system. 

[r(x) y' (x)] + P(x) y(x) = 0 (1. 1) 

y(a) = 0 y(b) = 0 

where r(x) and p(x) are positive , continuous, real valued functions 
on [a,b]. The system (1. 1) is known as a Sturm- LiouvilfeSystem. 

A frequent problem encountered in mathematical physics requires 
the determination of consecutive characteristic roots of a system (1. 1)l 
The solutions to (1. 1) usually do not possess the nice regular charac- 
teristics of the trigonometric functions; consequently the former are 
in general more difficult to manipulate. 

When the value of one characteristic root a of (1, 1) is known, 
Sturm's Comparison Theorem suggests a possible means of computing 
the next consecutive root b in the following sense: 

A sequence of functions is generated from (1, 1). Each element of 
the sequence is obtained by selecting a particular value of the indepen- 
dent variable x = x. and substituting this value into the coefficient P(x) 
and r(x) of (1.1) giving 

[r(xpy'(x)]' + P(x.) y (x) = 0 (1.2), 

The solution of (1. 2) is forced to vanish at y = a, and the next consec- 
utive zero of this solution is chosen as next element of the 

sequence. It is desired that the sequence so generated result in an 
equation of the form (1. 2) such that its solution vanishes at x = a 
and X = b, the zeros of the solution of (1, 1). 

The theory is developed for a general Sturm- LiouvilfeSystem but 
specific applications are limited to Bessel*s Equation of order n . 

(n a positive integer). 
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We shall define an indexing of the roots of a solution to a differen- 



tial equation as follows: 

Let rj be the first root > 0, and r^ the next root greater than 
r j. Furthermore, if r^ is the i th root ^ 0, then is the 

(i+l)st root > 0. 

If is the given root of a solution to a differential equation, we 

shall define r, , , as the next consecutive root of the solution, 
k+1 
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2. Mathematical Background, 



This section contains theorems and mathematical developments 

which will be called upon in later sections. 

2, 1. Sturm's Comparison Theorem. 

Given two equations of a Sturm- Liouville type; 

y» + Ay=0 (2.1.1) 

y(a) = 0 

u" + Bu = 0 (2. 1. 2) 

u(a) = 0 

where A and B are positive, real constants and B is greater than A. 

Let us determine a solution y(x) to (2. 1. 1) such that y(a) = 0. 

The elementary theory of differential equations tells us that the general 
solution is of the form 

y(x) = C j sin ^Ax + C 2 cos "^A x . 

When X = a 



y(a) = Cj sin”>y(^a + cos'Y^a = 0 
C j sin'^L a 



C2 = 



COS 






Since is an arbitrary constant, choose VTa. 

y(x) = [sin ~\fAx cos Va a - c OS VTx sin VTa] 

y(x) = sin (x - a). 

Similarly, (2. 1. 2) will have a solution u(x) = kj sin Vb(x - a) . 

The assumption A < B tells us u(x) will "oscillate" more rapidly than 

y(x) or that u(x) will have at least as many zeros as y(x) on a given 
• ..interval (a, b). 

When the coefficients of (2. 1. 1) are positive, real valued functions 

of the independent variable x, the solutions obtained must be defined 
in terms of so-called "special" functions such as Bessel functions. 
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Charles Sturm (1803-1855) made a study of the rates of "oscilla- 
tion*' of solutions to equation# of the form 
[ r (x) y'(x) ]' + P (x) y(x) = 0 
[ r (x) w'(x) ]' + Q(x) w(x) = 0 
y(a) = w(a) = y(b) = 0 

where r(x)_>0; P(x) and Q(x) are continuous on [a,b]. He proved 
that, in general, the larger P(x) the faster y(x) will oscillate. More 
concisely, Sturms Comparison Theorem is stated as follows: 

Theorem 2. 1 Sturm's Comparison Theorem. 

Given two equations; 

[ r(x) y']' + p(x) y = 0 (2.1.3) 

and [ r(x) v' ]' + q(x) v = 0 (2. 1. 4) 

where r(x) is positive in the closed interval [a,b]; 
r(x), p(x) and q(x) are continuous in [a,b] and 

p(x) with strict inequality for at least one point in 

[a,b]. If a solution y(x) of (2. 1. 3) has consecutive 

zeros at x = x and x = x, , with x, > x , and x, , x 
o 1 1 o 1 o 

are in the interval [a,b], then a solution v(x) of (2. 1.4), 
which vanishes at x = x^, will vanish again in the open 
interval (x^, Xj) . 

Proof: 

Suppose y(x) > 0 in (x^, Xj) and y'(x^)>0, y'(Xj) < 0 and 
v'(x^)>0. Then if y(x) and v(x) are solutions to (2. 1.3) and 
(2. 1.4) respectively, we have the identities 

[ r(x) y'(x) ]' + p(x) y(x) = 0 (2. 1. 5) 

and [ r(x) v'(x)]' + q(x) v(x) = 0 (2. 1.6) 
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If we multiply equation (2. 1. 5) by -v(x) and equation (2. 1.6) by 
y(x) we have 

-[r(x) y'(x) ]' V (x) — p(x) y(x) v(x) = 0 (2, 1. 7) 

[r(x) v'(x)]' y(x) + q(x) v(x) y(x) = 0 (2. 1.8) 

Adding (2. 1. 7) and (2. 1. 8) then integrating from x^ to x^ we 
have 



|Xi Xj 



•(x)[y<x)v'(x) - y'(x)v(x)] | +T [q(x)-p(x) ]y(x)v(x)dx = 0 

o J-x. 



or 



rx, 

r(Xj)y'(xj)v(xj) s / [q(x)-p(x)J y(x) v(x) dx (2. 1. 



9) 



o 

since y(x^) = 0 if x^ is a solution to (2. 1.3). 

Suppose now that v( x ) > 0 everywhere in (x^, Xj), (i. e. if 
v(x^) = 0 and v* (x^) > 0 then v(x) has not vanished anywhere in 
(x^, Xj) ), then this implies that r(xj) v(x^) y'(xj) is negative while 



X 

(q(x) • p(x) ) y(x) v(x) dx is positive. This contradiction implies 
o 

that v(x) cannot be positive everywhere in (x^, x^). Similarly> the 
assumption v'(x^) < 0 and v(x) negative everywhere on (a,b) again 
leads to a contradiction and the theorem is established. 

2.2 The Normal Form. 

Consider a differential equation 

y" + P(x) y' + Q(x) y= 0 (2. 2. 1). 

Let y = u V, then y* s uv' +u'v 




and 



y‘% uv'* + 2u'v' + u"v. 

Substituting lor y in (2.2.1) we have 

uv" + [2u' + P(x)u]v* + [u" + P(x) u' + Q(x)u] v s 0 (2.2.2). 
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Suppose we choose u such that the coefficient of v* vanishes. Set 



2u* + P(x) u = 0 
du/u + dx = 0 

u = C exp [- 

Let C = 1 and we have 
u* = ^ P{x) u 

u” = ^P(x) u* - ^ u ^^^dx 

^ (P(x))^ 1 dP/j 

= + — - / dx u . 

Substituting u' and u" in (2. 2. 2) we get 



uv 



" + [-P(x) u + P(x) u] v' + [Q(x)u + u - i ^^/dx u 



P(x)"' 1 . 

|-i- u J V = 0 . 

We see that we have a factor of u 

divided out leaving 



in every term, therefore it may be 



v" + [Q(x) - ^ ^ ] V = 0 (2. 2. 3) 

We say that (2. 2. 3) is equation (2. 2. 1) transformed into normal 
form. We will see that the normal form plays an important role in our 
later development. 

2. 3. Existence Theorem. 

Consider a differential equation reduced to normal form 



y" + R(t)y=0. (2.3.1) 

y(a) = y(b) = 0. 

We desire to know if there exists a t in the interval (a,b) such 
that the equation 

u" + R(F) y = 0 (2. 3.2) 

R(F) is R(t) evaluated at t = F such that a solution u{t) to (2, 3,2) 
can be forced to have zeros at both t = a and t = b. The answer to 
this question may be stated in the form of a theorem as follows: 
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Theorem 2. 2. 



Given a differential equation of the form: 

+ R(t) y = 0 (2. 3.3) 

where R(t) is positive, real valued, continuous and monotone in (a,c), 
where b < c, and whose solution y(t) has consecutive zeros at t = a 
and t = b, a < b < c. Then there exists a t in the interval (a,b) such 
that the equation 

u*» + R(F) u = 0 (2. 3.4) 

has a solution u(t) with a zero at t = a and moreover, u(t) has its 
next consecutive zero at t = b. 

Proof: 

First assume R(t) is monotone decreasing, then 

R(a) ^ R(t) > R(b) for all t in (a, b). Choose that solution u^(t) of 

u** + R(b) u = 0 such that u^(t) vanishes at t = a. A solution that 

satisfies this requirement is ^(t) = sin y^R(b) (t-a). Clearly 

u^(a) = 0 and moreover, u^(t) will have its next consecutive zero, 

for t greater than a, when the argument of the sine function is equal tt, 

that is , when VR(b) (t~a) = ir . Since R is assumed decreasing, 

R(b) < R(t) for all t in (a, b[ therefore, Sturm*s Comparison Theorem 

tells us that will vanish again at t = t^ where b < t j. Similarly, 

let u** + R(a) u = 0 have a solution u (x) which vanishes for t = a. 

a 

Then with R(t) < R(a) for all t in the interval (a,b), have that 
u^(t) will vanish again at some point t = t^ such that a < t^ < b. 

Case II. 

Now suppose R(t) is increasing. Choose u^(t) a solution of 
v" + R(b) V = 0 

such that a root at t = a. R(t) < R(b) for all t in (a,b), hence v^(t) 
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will vanish again when t = < b. Similarly, v" + R(a) v = 0 has a 

solution V (t) which vanishes at t = a and will have its next consecutive 

ci 

zero at some t = t., b < t., 

4 4 

In either case, R(t) increasing or R(t) decreasing, the root b of 
equation (2. 3. 3) is between two roots of (2. 3. 4) determined by evalua- 
ting R(t) at t = a and at t = b. This and the continuity of R(t) assures 
us that there exists a t in (a,b) such that ut^ (t), a solution to (2. 3.4) 
which vanishes at t = a, will vanish again when t = b, and the theorem 
is established. 

An interesting conclusion from Theorem 2. 2 is the following: 
Theorem 2.2. 1. 

Given a system of the type (2. 3. 1) where R(t) is monotone, then 
the next consecutive root of a solution u(t) to (2, 3.4) such that 
u(a) = 0 must lie in the interval (a + , a = ) . 

y R(b) Y^) 

The proof follows directly from the proof of Theorem 2. 2. The end 

points of the interval (a + ^/yR(b), a+ >/r(I) ) are precisely the 

t. found in Theorem 2. 2. 

1 

Note that the arbitrary constants occuring in the solutions to equations 
(2. 3. 1) and (2. 3.4) have no effect upon the zeros and, therefore, will 
be ignored in all further considerations. 

We should also note that in a general Sturm- Li ouvilh System (2. 3. 1), 
the function R(t) is not always monotone. 
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y”+xy* + (x^- n^) y = 0 



3. Bessel’s Equation. 

Let us now consider a typical Sturm -Liouville System: 

(3. 1) 

y(a) r: y(b) = 0 (3. la) 

Equation (3. 1) is the general form of Bessel’s Differential Equation of 

order n. A solution of (3. 1) is 

y(x) = J^(x) (3. 2) 

where J (x) is known as the Bessel function of the first kind, 
n 

We see that before we can derive any benefit from the theorems of 
section 2, we must reduce (3. 1) to its normal form. Write (3. 1) in 
the form 

(3. 3) 



y" + iy’ + (i-^)y = o. 

X 



Section 2. 2 tells us that (3. 3) may be reduced to normal form by 
the transformation 

y = V exp [ - j dx ] 

= V exp [ - 2 logg X ] 

(3,4) 



= V X 



-1 /2 



Then 



, 1 , -3/2 . , -1/2 

y' = - y v(x) + V X 



(3.5) 



y» = ^ V X 



-5/2 1 , -3/2 1 , -3/2 „-l/2 

: - v’x - v'x + v”x 



(3.6) 



Substituting (3.4), (3.5), and (3. 6) into (3.3) and simplifying we 



have the desired normal form of Bessel's Equation 

l-4n^ 

v" + ( 1 + — ^ ) V = 0. 

4x 



(3. 7) 



The transformation (3. 4) implies that (3, 7) has a solution 
v(x) = x^ J„(^)* (^* 



n 



Trivially any zero of (3, 2) is also a zero of (3„ 8). 
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Consider the differential equation 

+ w = 0 (3. 9) 

w(a) = 0 

We know that a solution to (3, 9) 
w(x) = k sin(x-a) 

vanishes at x = a and again at x = a + ir . 

Suppose we examine (3.7) and (3.9) with respect to Theorem 2. 1 
(Sturm' s Comparison Theorem). Let us denote the coefficient of v in 
(3. 7) as 

P(x) = 1 + { 3 . 10 ) 

4x 

We see from (3. 10) that we must consider two cases, 1) when n = 0 
and 2) when n 1 . 

Case I. 

n = 0. We see that P(x) > 1 and continuous, therefore, theorem 
2. 1 applies. We conclude that if the solution (3. 8) to (3. 7) vanishes 
when X = a, it will vanish again on the interval (a, a + it). This implies 
that the difference between consecutive roots of J^(x:) is less than ir . 
Moreover, we see that 
lim P(x) = 1 . 

X 0=5 

This tells us that with respect to the index i of roots r^ ; 
lim 1 I = IT . (3. 11) 

Case II. 

n 1. In this case 0 < P(x) < 1 when 7 4n -1 <x< «e^ . 

Then x must be in the interval ( ^ V4n^ -1 , o® ) in order that (3, 7) 
have an oscillating solution. Theorem 2. 1 tells us that a solution to 
(3. 7) which vanishes at x = a will not vanish again in (a, a + it ). 
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Therefore, the difference between the consecutive roots of J (x) is 

n' ^ 

greater than *rr . Similarly, 

lim P(x) = 1, when n ^ 1 , 

tells us that the limiting difference is the same as in Case I, namely the 
difference given in (3, 11)« 

The treatment of the differences between consecutive roots of the 
Bessel functions in the above two cases, suggests that further treatment 
of the problem also be made in two parts, i« e„ 1) n = 0 and 2) n ^ 1 
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4. Application to Bessel's Equation of Order Zero« 

When n = 0, the nornaal form of Bessel's differential equation 
reduces to 

v" + (1 + — K- ) V = 0, (4. 1) 

4x 

Suppose the conditions 

v(a) = v(b) = 0 (4. 2) 

are imposed. Let the coefficient of v in (4. 1) be denoted by P(x), i. e, 

P(x) = 1 + (4. 3) 

4x 

Clearly P(x) 1 for all x, P(x) is continuous and decreasing^ Then, 

from theorem 2. 2, there exists an x in the interval (a,b) such that 

the equation 

u" + (1 + — ^ ) u = 0 

4x 

has a solution 

u (x) k sin 

with a root at x a and, moreover, will have the next consecutive 
root at X = b. 

Let us now attempt to construct a sequence which will converge 

to X « There are two immediate choices for the first element of the 

r I' 

sequence, a or a + tt « We will generate the sequence by means 

of the following recursion formulae: 



P(x ) - 1 + 


(4.4) 


4x. 




1 




u!* f P(x.) u. = 0 

1 ' i' 1 


(4.5) 


u.(x) = k sin V^P(x^)(x-a) 


(4.6) 


^i+1 ~ ^ 


(4. 7) 
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Wherein (4, 6) is a solution to (4. 5) which is forced to have a root at 
X = a; and is the next consecutive root of (4.6). We see from 

(4.4) that P(x^) ^ 1 for all x. . This implies that < a + ir for 

all x^ . 

Case I. 



Choose x^ = a, giving 



P(Xq) - 1 + — = 1 + ^ 

4x 4a 

o 


(4.8) 


u‘* + P(x ) u =0 
o o o 


(4.9) 


Uo(x)*= sinyp(x^)(x-a) 


(4. 10) 


Xj = a + ir /y"P(x^) 


(4. 11) 


Since P(a) P(x) for all x in the interval 


(a,b), theorem (2.1) 


(Sturm's Comparison Theorem) tells us that 


Xj lies in the interval 


a < Xj < b , 




Then Xj is substituted in the recursion equations giving 


p(Xi) = 1 t ^ 
4xi 


(4. 12) 


u'l + P(xj^) Uj = 0 


(4. 13) 


Uj^(x) - kj sin ^P(xj^) (x-a) 


(4. 14) 


x^ = a + IT / yP(x^) 


(4. 15) 



We cannot apply theorem 2. 1 with respect to (4. 13) and (4.3) 
since we cannot guarantee the necessary inequality between P(Xj) and 
P(x). However, we may apply theorem 2.1 with respect to (4. 13) 
and (4. 9) since P(xj) < P(x^). This tells us that x^ ^ Xj . 
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If we expand \4o 15) 



X 2 = a 





we see the sequence is generating a rather strange appearing continued 
fraction; strange in the sense that the more common continued fractions 
do not contain radical expressions* Moreover, note that the fractional 
form becomes exceedingly more complicated as each new term is added* 
Let us now assume that for some m >0, we have 

*m-l > ’‘m-2 

The inequality (4„ 16) implies that the following inequality also holds 

(4. 17) 



P(x J = 1 + -V- < 1 

m- 1 . ^ m-^: 4 Y ^ 

m-1 m-^ 



Theorem 2* 1 immediately tells us that x is greater than x . 

^ m ® m-1 



Inductively we have shown that the sequence is increasing* 
Moreover, we have shown that the sequence is bounded above by a + tt, 
consequently, the sequence converges to some ic such that 

TT 



X == a + 



1 + 



1 

4x 



t'Tt? »2 ^ c'/’ 2 2 

(x-a) - 4x TT 



4x > 1 



(4x 4 l)(x - a)^ 



^ 2 2 y-w 

4x TT = 0 



4x ^ - 8a X ^ + (4a^ + 1 - 4 it^) x ^-*2a x + a^^= 0 * 

Then x is the root of (4. 24) that is nearest a + it 

^ 2 , 



(4„23) 

(4o 24) 

We see from 



(4, 23) that for large roots, 4 x 1 and the large roots of (4* 24) may 
be approximated by: 
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4 *^'" ^ ,'i^ 2: 2a - 

X (x-a)-4x ir=0 

/•- x2 2 n 

(x - a) - TT =0 

O 

X = a + TT 

A simple calculation will show for any given root *a* of J^(x) that 
(4. Z4) has a root near a + tt . 

Case IL 

Let us choose x = a 4 tt , 
o 



P(x ) = 1 + 
o 



1 



4(a+ irf 



u*‘ + P(x ) u =0 
o ' o' o 



u (x) = k sin '/pix T (x-a) 



(4.25) 

(4.26) 



Xj = a + TT / 






(4.2 7) 

(4. 2 8) 

We have P(a + tt) < P(x) for all x in (a^b), therefore, we know x 
cannot be in (a,b), therefore, we see 



1 



b < Xj < a Hh tr 

Substituting in the recursion formulae we have 

1 



P(Xj) = 1 + 



4x 



1 



Uj" + P(Xj) 



Uj(x) = kj sin '/p(x^) (x-a) 
x^ = a + -rr / V^P(x^) 



(4. 29) 

(4. 30) 

(4. 31) 
(4.32) 
(4. 33) 



The inequality ^4«29) tells us that P(xj) ^ P(x^), and we get < 'k.j 
from Theorem 2„ L Moreover, (4« 29) also tells us that P(xj^) < P(x) 

for all X in (a,b) and ^2 ^ • 

If we reverse the assiimed inequality in {4« 16) we can show induc- 
tively that this sequence £xj'' , generated by starting at - a + tt , 
is a monotone decreasing sequence^ Moreover , we see that |x^| > b 
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for all i, hence is bounded below. Thee the sequence converges to an 
X|^ where 

Xj = a + TT / '/p(x^ ) 



tr 




(4 + l)(xj-a)^ - 4 TT^ - 0 

4 - 8a Xj^ + (4a^ + 1 - 4tr^) - 2a - 0 (4« 34) « 

The quartic (4« 34) lias coefficients which .are identical to the coef- 
ficients of like powers of in (4,24) ^ therefore we conclude that 
Xj = X , This shows that there is no difference in the limit between 
the sequence generated in Case I and that generated in Case II, Un- 
fortunately, we have also shown that x > b , but we have no indica.tjOn 
of the difference between x and b. 

The quartic (4, 24) has been solved with the aid of the CDC 1604 
Digital Computer for the first 20 roots of J^(x), The computed values 
of X are compared with the tabulated values of the roots of « The 

difference between x and the tabulated value of b is shown in the error 

column. The difference has its maximum at the first root and decreases 
quite rapidly as x grows larger. 

Simple calculations show that ^ is almost midway between b and 

a + u for the first root, but this relation does not persist as the root 
number increases. In fa.ct x approaches b more rapidly than a + tt , 

The fact that the sequence x, does not converge to b implies that 

some change to the recursion formulae must be made. The most direct 

method of changing the formulae consists of a modification to P(>:)„ 
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4. 1. Linear Average. 

Suppose a sequence were to converge such that (4. 1. 1) holds, 
b = a + ir / Y 1 + — 

9 A 



1. e. 



(4. 1. 1) 



4S 



where %• is the limit point of a sequence 



' ^ + l)(b-a)^ = A% 

4% ^[(b-a)^ - ir^] = - (b-a)^ 



* = ^ 



(b-a)^ - iT^ 



% 



b-a 



vt 



^ (b-a)^ 



We have shown that (b-a) < tt , therefore, we know that % is real« 
Computed values of % are displayed in table II, along with the tabulated 
values of b. The error column is the difference between x and M , 



and indicates the amount x must be decreased if the sequence is to 
converge to b. Noting that % is almost midway between a and x the 

following modification to the recursion formulae is suggested: 

1 



P(x.) 


= 1 + - 2 
4(a + X.) 


(4. 1.2) 




— Z 




u. f 
1 


P(x.) u. = 0 
1 1 


(4. 1.3) 


u.(x) 


= sinV^(x-) (x-a) 


(4. 1.4) 


^i+1 


s a IT / 


(4. 1.5) 



In essence, the modification consists of merely substituting 
{a X ) ^ 

' i' in place of x. in {4«4)« Clearly the conclusions reached 
— " 

regarding convergence of the sequence | are also valid for the 
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sequence | ^ • For the sequence y we are unable to make any con- 
clusion relative to b since the of Ca,se II is 



X 



1 



a + TT / 





a 



+ 



TT 




"7 

) 




(4. 1.6) 



and lies in the interval (a,b). Moreover, we can see that for i ^ 1 
each element of the sequence is smaller than the corresponding 



element of the sequence | we conclude that the limit of the se- 



quence 






X is smaller than x . We have that 



X = a + 






= a + 



. ,a+x» 
4 {^— ) 




Vl + 



(a + x) 



(x - a)^[ (a + x)^ + i ] - (a + x)^ = 0 (4. 1. 7) 

X + (1 - 2a^ - X [2a(l + ir^)] x + + a^(l-ir^) = 0 (4. 1. 8) 

For the larger roots of (4. 1. 8) we look at (4,, 1. 7) and see that 

2 

(a + x) 1 so that an approximation yields 

(x - a)^ (a + x)^ - (a + x)^ ir^ = 0 
X = a + ir . 

We conclude that (4. 1. 8) has a solution near a + tt, and that this 
solution is the x we desireo 

Table III lists the desired roots of (4„lo8) for the same 20 roots 
as in section 4. A comparison of the error column of table III against 
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table I shows that the difference between x and the tabulated root b is 






about one-tenth as large as the difference between x and b at the first 






root. Moreover, we see that x approaches b more rapidly than x, 
that X - b is only . 01 times ^ - b for the tenth root. 

B ased upon the results of these computations, we conclude that 
although X is nearer b than x , x is still not the value of x whose 
existence is guaranteed by theorem 2. 2* 

4.2. Mean Value. 

Consider the integrated mean value of P(x) over the interval (a,x) 
and denote this integrated mean value by P(x), then 

^ X , 



P(x) = 



= — f 

Ja 



(1 + 



4t 



) dt 



1 



x-a 

J 

x-a 



Ft - — 1^ 

^ 4t ■'a 

[x-a + T ( - - r ) J 



1 



= 41 ^ ) 



= 1 + 



4ax 



We see that 



(a+x)'^ 

i V 1 

■ (a+x)^ 

a^ + 2ax + x^ ^ 4ax 
a^ - 2ax + x^ ^ 0 
(a-x)^ > 0 for all x 4 

Suppose we generate a sequence 
recursion formulae: 



(4.2. 1) 



a . 

by means of the following 
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1P(X.) = 1 + ■— 

' 1 4ax* 

1 


(4. 2.2) 


u.** + IP(x.) u. = 0 


(4.2.3) 


Uj(x) = k. sin 'Y^(x-)(x-a) 


(4.2.4) 


^i+1 = a + ir/’\/BP(x.) 


(4. 2.5) 



But equations (4*. 2„ 2) through (4, 2. 5) are the same as the set of 

2 

equations (4.4) through (4.7) except ax^ is substituted for x^ in 
(4.4). If we make this substitution throughout our preceeding considera- 
tion we see that the sequence also converges to some x = x such 

that 




(4ax + l)(x-a)^ - 4ax tt^ = 0 

4ax^+ (l-8a^) + (4a^ - 2a - 4a tt^) x + = 0 (4. 2.9) 

Inequality (4. 2. 1) shows us that x must indeed be less than x 
but once more we are unable to state any analytical conclusions regard- 
ing X and b . 

Table IV contains the roots of (4. 2. 9) for the same zeros of 
J^(x) as previously computed. We see that x is smaller than b. We 
note also that |x-bj is greater than |x-bj , 

Since b has been bracketed by x and x , a linear average of the 
two is suggested. Accordingly we define 

P(xl.) ^ (P(x.) + P(x^)/2 (4.2.10) 

xL^ j - a + -it/ Vp(xI^) (4. 2. 11) 

Then equation (4. 2. 10) and (4. 2. 11) together with equations (4. 5) and 
(4. 6) constitute a set of recursion formulae which will generate a 
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sequence |xl.' ^ We know the sequence xi will converge to some 

r ^ .j. > 

xl since it is the sum of two convergent sequences arid ^ x/,' . The 

values of xl are displayed in Table We see that although the values 

of xl are nearer to b than they are not as near as x « 

o— 

Suppose we consider the average of x and x . Let us define: 

P{x2.) = (P(x.) + P(x.))/2 (4.2.12) 

^^i+1 “ a + it/ V"p(x2^ ) (4.2.13) 

Again equations (4.2.12) and (4,2.13) together with (4.5) and (4.6) 

will generate the convergent sequence /x2.(f^ . We see from Table VI 

that x2 is nearer to b than any of the previous values generated^ The 

“ -3 

maximum difference between x2 and b is less than 10 for the first 

7 

root and decreases to less than 10 for the 20th root. The difference 
decreases quite ra.pidly as is evidenced by a value less than 10"^ beyond 
the fifth root. 

4.3 Cumula.tive Effect. 

Let us now investigate the cumulative error for each of the sequen- 
ces previously discussed. The errors in this case are cumulative in 
the following sense: 

Assume that a = r^, the first root of J^(x) is known, we will 
compute an x from a sequence. Then each successive root is generated 
using the root last computed as the value of the known root. 

The cumulative errors for x, x , x, xl and x2 are shown in 
Tables VII through Xl respectively. As well expected, the cumulative 
error is due mostly to the error in computing the first root. We see from 
Table XI tha,t the cumulative error for the 20th root is only 10% 
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greater than the error for the first rool^ The worst case is shown in 
Table X where the cummulative error ior the 20th root is 40% greater 
than the error for the first root« 
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(5. 1) 



5o Application to BesseVs Equation of Order n ^ 1 . 

When n J>. 1, the normal form of Bessel's equation is 

y»t + (1 + ) y = 0 

4x 

Let us impose the condition 
y(a) = y(b) = 0 . 

Again let P(x) denote the coefficient of y in (5. 1), i.e, 

1 -4n^ 

P(x) = 1 + (5.2) 

4x 

We see that P(x) > 0 for all x moreover, P(x) is 

monotone increasing, lim P(x) = 1, i*e« 0 < P(x) < 1 for all x 

X j , — - 

such that x is in the interval y4n^-l, . We shall limit our dis- 

cussion to the case P(x) > 0 in order that the conditions of the theorems 
of section 2 may be satisfied^, 

Except for P(x.), the recursion formulae defining the sequence will 
be the same as in section 4, i,e. 



P(x ) = 1 + 


(5. 3) 


4x 




o 




u!' + P(x.) u. = 0 

1 ' i' 1 


(5.4) 


u^(x) = k. sin^V^(x^) (x-a) 


(5.5) 


X. . , = a + ^/ Vp(x.) 

1+1 ^ ' i' 


(5.6) 



As before we consider the two cases x = a and = a + ir , 

o o 

Case I: 

Choose X = a , giving 



Ti/ ^ _ 1 . l-4n l~4n 

^ , Z“ - ^ + ~T 

4x 4a 

o 


(5. 7) 


u*' + P(x ) u =0 
o ' o' o 


(5.8) 


u^(x) =• sinVP(x^) (x-a) 


(5.9) 


Xj - a + IT / \/ P(x^) 


(5. 10) 
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Since P(x) is increasing, we see that P(x^) < P(x) for all x in 

(a,b). Then theorem 2. 1 tells us that b . Moreover, 1 > 

therefore, from (5, 10) we see that > a + tt , 

When X, is substituted in the recursion equation we have 
^ 2 

P(x ) = 1+ ^-l4- (5.11) 

4Xj 

Uj* + P(xj) Uj = 0 

u^(x) = kj sin Vp( X ) (x-a) 

= a + ir / P(xj^) (5. 12) 

We know from the inequality x^^b that x^ is in the interval (a,b). 
Trivially then x^ < Xj , moreover, P(Xj) < 1 tells us that x^ > a + tt « 
Then we have 



X =a<a+7r<x>^<b<x, 
o 2 1 



(5. 13) 



Let us now consider the next element of the sequence. If x^ is 



substituted in equations (5« 2) through (5, 6) ,we have 
P(x ) = 1 + 



(5. 13) 



u 



2 + ^<^2^ ^^2 = ° 



u,(x) k, sin Vp(x,) (x-a) 



Xj = a + tr/Vpiix^) 



‘2'''' ■ ■'2 ° ' 2' 

(5. 14) 

Since is in the interval (a,b) we do not have the necessary 

inequality for theorem 2, 1 and consequently we can no longer determine 
the relation of x^ to b for i ^ 3 , We see, however, that PCx^ )>P(Xj) 
since P(x) is increasing and x^ < Xj , This implies that x^ ^ x^ • 
Clearly PCx^) < P(x^) and < Xj . Then we have 



X = a < a + IT < < X., < X| 

o ^,.6 1 



(5. 15) 
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Suppose, for some ^^3 that the following inequalities hold: 



a < a f « o < x, ^ < x-, < x, . < 

^ k- Z k k - 1 



< X 



The recur sion formula gives 



IT 
k-2 
2 



P(x^) 



1 + 



P(Xk_i) = i + 



4 X 
J,-4ri 

I -4n 



(5, 16) . 

(5. 17) 
(5. 18) 
(5, 19) 



4x. 



k-l 



^^k-2 ^ ^^^k.2)"k=.2 " 0 


(5. 20) 


u|^ + ° 


(5.21) 


^k-1 ^^^k-i^^k-i “ ° 


(5.22) 


1 = a f tr / Vp(Xk 2 ) (x-a) 


(5.23) 


Xk = a + TT / Yp{Xk^^) (x-a) 


(5. 24) 


Xk^j = a f ir ri P(Xk) (x-a) 


(5. 25) 



The assumption (5„ 16) tells v- s that Pix^ 2 ^ ^^tXk)^P(Xk ^) since 
P(x) is increasing. If we apply theorem 2, 1, first with respect to (5. 21) 
and (5,20), then with respect to (5, 2i) and (5„ 22); we have 

= a < a f w < >2 < o . , < x^ 2 ^ ^k+l*" 1^" ” ° *^^1 

We have now shown inductively that lies in the interval 



(x. , X. ,) for all 1 . Let us now select the sub sequence |z .1 of |x. ,? 

11-1 

such that for each z. - « Trivially “x „ Assuming 

j 1 00 

for some m^O we see that z 



then z y z 
m 



m 






^ m- 1 

Moreover^, we knew that > lies in the interval (x, , x., 

Ci mir 1 ; tj m Ca m i i 



therefore a we conclude z 



mf"! 



y ^ 



m 



and that if z .!> is monotone increasing^ 
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Clearly, the sequexice is bounded abo^/?i Dy Therefore, we 

conclude that the sequence Zj converges ro some z „ 



Similarly, consider the subsequence of chosen such that 

T' 1 

for each n, argument similar to the one for 'we 

can show that the subsequence is decreasing and bounded below by 

a, therefore, it converges to some v « We have now shown that the 



sequence I x/| contains two convergent subsequences JZj ^^and ;v^ ■ , and 
in order to prove convergence of we must show^ that the subsequences 
converge to the same point, i« Co that z ^ We know from (5„ 25) 
that 



a 4- 






a ¥ 






w 


C 








4 /a + 



(5„ 26) 



ir 

/iT il! 

I 4 F ‘ 



and V 



a Ht 



/ 



I -in 
4 / a t 



{X 27) 



\ 



1 f 



1 -4a ^ 

7T^2 I 
4v / 



Equation (5,26), after considerable algebraic manipulation, can be 
reduced to the equa.'Ion 

[z-a)^^(4 I* a(4z " f l-4n‘ ) f 4z ""w’" | + f I-4n ' )(4z 1 -4n^')«4ir'^ 

;^a(4z -f l«4n ) 4z tr | ] 1 6a'"z^(4z*'‘ir l-^'4n'' j , (z^ a)~- 4Tf -» 0« 

which is clearly ai polynoimal of degree 8 „ Similarly, Equation (5. 27) 
will reduce to the same poiynomiaJ, Then we can conclude tha,t z v 



and the sequence Iconverg^-- s to some x such that 
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^ I 

X == a + 



which reduces to 



TT 



1 + 



1- 4n 



4x 



4'x ^-8ax ^ + (4a^+ l«4n^-4ir^)x ^-2a(l-4n^) x 

+ (l-4n^) = 0 (5. 29) . 

Moreover, x is the root of equation (5« 29) nearest a + tt . 

Case II: 

Suppose we choose = a + ir then (5„ 3) through (5. 6) becomes 

P(x ) = 1 + = 1 + ^ (5. 28) 

^ 4x^ 4(a+7r) 

u** + P(x ) u - 0 
o o o 



u^(x) = sinV P(x^) (x-a) 

X^ = a + TT / V 



( 5 ^ 29 ) 



We have chosen an x^ with interval (a,b) , therefore, we cannot 

make any comparisons between x^ and b. Since 

P{x ) < 1 we know x, > a + tt = x 
' o 1 ' o 

Let us now find the next element X 2 of the sequence . We 

have the equations 



p(-i) - ^ 



2 



4x 



1 



u'^' + P(vj) Uj = 0 
Uj^(x) = k j sin Vp(xj) (x-a) 




s therefore, tells us that < Xj^ , 



The inequality x, ^ implies that P{X j) < P(x^) since P(x) is 
decreasing. Theorem 2. 

moreover, from P(> ,| < 1 we see that x^ ^ a, + tt giving 

(5.30) 
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X 



a + 1 T < X 2 <• X j 



•ssBies- 



I 



m 








Clearly* ior ail i> Z this sequence hv-^Laves exactly as the sequence 



of Case I, therefore* converges to some Xj such that 

/■ji? 

“ a + 



1T 



1 



1 + 



1 -4n 

. ^ 2 
4x 



(5.31) 



But equation (5»3^1) after simplification reduces to the same quartic as 
(5o 2 7), and we conclude that = x « 

In either case for a > 1 we have been unable to relate x to b, 
therefore, we must consider computations for the first 20 rootso 

The results of the computation displayed in table XII are com- 
parable to those shown in table T for the case n s 0, as should well be 

expected. We see that for the 20th root, the error in table XII is 

-4 - 6 

10 while for J , it is about 5 x 10 
o 

Suppose, however, the sequence to converge to b. The 

final result would for some %, the following 

b a f tr 

i/ptn 



= a h 






1 4 



U4n 

4^^ 



(4X ^ t i. 4o“) {b-a}- - 4:i ^ tr^ 

which when solveo lor % snows that 




ir^) 



(5. 32) 



Trivially t. d*4»irtnJafed in this manner is real »-aiued, since 
(b-a)^-tr and 4n^ ^ 1 . The values of % computed from equation (5.32) 
are shown in table XIII. These values suggest that the following modifi- 
cation be made to the recursion formulae. 
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m 



• i 




•> 



P(x.) 



i + 




u " + P(x.)a = 0 

1 13 

u^(x) - sin VlP(x^){x-a) 

^i+1 “ ^ 

We have already shown that the sequence ,^x. f defined by equations 
(5.3), (5.4), (5„ 5) and |5„ 6) converges to x. The sequence is 

the sum of a constant and a multiple of x^, and therefore will converge. 
This implies that the sequ'^nce | |converges to some x such that 




x = a + 



\ + 



4(2$)' 



(6.13) 



x^+(l- 4 n^- 2 a^-Tr^")x" - [ 2 a|i * 4 n^+ir^) ]x 4 i.*-i*i* - ir^) - 0 ( 5 . 14 ) 

The roots of the quartc (So 14 ) nearest a r ?r are displayed in 
table XIV for th< first roots of Jj(x)o Comparing the errors of 
table XIV against U<ble XU we see that jx>bj is one order of magni- 
tude smaller than j^^ bj o 

We shall 'amine the sequence |x^ | i, io^o the mean value of 

P(x) over the interv,#! (a^ The mean value oi P(>) is defined in the 



same manner as in clion 4^ 2^ We know that the sequence | x. I will 



l 1 , 



converge since it consists of elements which are multiples of the elements 
of a convergent sequent ^ o Then we have 



= a + it/ y P(x 

ir 



a 4 






1 + 



1^4n 
4a X 
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A 



(5, 15) 



4a^^+(l-4n^- 8a^)x'^+(4a^-2a«8a^"n^-4aTr)x a - 0 

We see that for n ^ 1, the limit of the sequence results in a cubic 
similar to the one for x when n = 0, Again we know that x is the root 
of the cubic nearest a + tt. 

Table XV displays the roots of equation (5„ 15) for the first 20 
roots of Jj^(x). As expected from the results for Jq(x), we see that 
again x is past b and that x and x have straddled b. We note, 
as for Jq(x), Ix-bj > jx-b] . 

We shall now examine xl and x2, as defined in section 4. 2. The 

values of xl are displayed in table XVI , Again, as for J^(x), we see 

that jTT is further from b than x. 

is again nearest b for Jj(x), with a maximum error for the 

first root. In general, we see that the vaiious sequences generated for 

Jj(x) tend to converge to points farther from b than the same sequence 

for J (x). 
o 

As expected, the cumulative errors shown in tables XVIII through 
XXII exhibit the same characteristics as the cumulative errors for 

Jo(x). 
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6. Conclusions. 



We have proposed a technique of computing the consecutive zeros 
of a solution to a Sturm- Liouville System. 

[r(x)y'(x)]' + P(x)y(x) = 0 (6. 1) 

y(a) = 0 y(b) = 0 

The technique requires the determination of an ^ such that the equation 

[r(x)u'(x)]' + P(x)y(x) = 0 (6.2) 

will have a solution u(x) which maybe forced to vanish at x = a and 

at X = b . 

The existence of such an x has been shown by theorem 2. 2. More- 
over, in theorem 2. 2« 1 we have established upper and lower bounds on 
the range of this x . 

In the specific case of the Bessel functions* we note that each of the 

sequences generated by the proposed techniqu-^ co ./erges to a particular 

X, no matter whether the first element was chosen equal to the known 

root a, or chosen equal to a + tt . We see that x turns out to be the 

root of a quartic nearest a + tt. Similarly^ we see tha.t x is also a root 

/\ 

of a quartic polynomial while x resolves into a cubic equation. 

The solution for % for the case n = 0 and for n = 1 both indicate 
that the value of rht; desired x as defined by theorem 2^ 2 is almost 
midway between a and bo The difference between x and % approaches 
tr/2 as x incr }ases. 

We have also shov/n that the desired x for the Bessel functions lies 
between the x determined by a linear averaging technique and the x 
determined by the mean value of P(x) over the interval (a,x)o 

Unfortunately, all of the above generating techniques failed to result 
in a sequence converging to the desired x. 
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TABLE 1 



VALUES OF X *** Jq(X) 



ZEROS 


OF THE BESSEL FUNCTION OF ORDER ZERO 


FOR ROOTS 


ONE THROUGH TWENTY, BASED ON THE RECURSION 


FORMULAE 


P(X(in = 1 + 1/4(X(I > 


X(I+1> » 


A+TT/ \/p(x( in 


ERROR 


IS THE DIFFERENCE 


BETWEEN TABULATED 


AND COMPUTED ROOTS 


ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


5.5200781103 


5.5336719233 


.0135938131 


2 


8.6537279126 


8.6564432415 


.0027153292 


3 


1.1 .7915344390 


1 1.7925004759 


.0009660370 


4 


14.9309177082 


14.9313671612 


.0004494532 


5 


18.0710639674 


18.0713085639 


.0002445965 


6 


21.2116366294 


21.21 17842031 


.0001475737 


7 


24.3524715304 


24.3525673207 


.0000957905 


8 


27.4934791317 


27.4935447974' 


.0000656658 


9 


30.6346064676 


30.6346534281 


.0000469606 


10 


33.7758202134 


33.7758549489 


.0000347356 


n 


36.9170983527 


36.9171247659 


.0000264135 


12 


40.0584257636 


40.0584463133 


.0000205500 


13 


43.1997917127 


43. 1998080136 


.0000163012 


14 


46.3411883712 


46.3412015187 ' 


.0000131483 


15 


49.4826098960 


49.4826206546 


.0000107595 


16 


52.6240518400 


52. 6240607537 


.0000089143 


17 


55.7655107528 


55.7655182229 


.0000074702 


18 


58.9069839241 


58.9069902431 


.0000063197 


19 


62.0484691886 


62.0484745828 


.0000053946 


20 


65. 1899647992 


65. 1899694391 


.0000046415 
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table II 

VALUES OF **• Jq(X) 



ZEROS 


OF THE BESSEL FUNCTION OF ORDER ZERO 


FOR ROOTS 


ONE THROUGH TWENTY, BASED 


ON THE RECURSION 


FORMULAE 


CHI 


IS THE VALUE X MUST 


ASSUME TO CONVERGE TO 8 


50 = ( (8- 


•A)/2)( Vi/((B-A)^< 


-7i2, ) 




ERROR 


IS THE DIFFERENCE 


BETWEEN X AND CHI 




ROOT 

NO, 




TABULATED 

ROOT 


COMPUTED 

ROOT 


ERROR 


1 




5.5200781103 


3.8368678503 


1.6968040731 


2 




8.6537279126 


7.0180586502 


1.6383845913 


3 




11.7915344390 


10. 1751194785 


1.6173809974 


4 




14.9309177082 


13. 3249362609 


1.6064309003 


5 




18.0710639674 


16.4716300704 


1.5996784936 


6 




21.2116366294 


19,6166964089 


1,3950877941 


7 




24.3524715304 


22.7608070695 


1.5917602512 


8 




27.4934791317 


25.9043037482 


1.5892410493 


9 




30.6346064676 


29. 0473744823 


1.5872789458 


10 




33.7758202134 


32. 1902150838 


1.5856398651 


11 




36.9170983527 


35.3327245601 


1.5844002059 


12 




40.0584257636 


38.4751889762 


1.5832573371 


. 13 




43.1997917127 


41.6175294435 


1.5822785702 


14 




46.3411833712 


44.7596824104 


1.5815191083 


15 




49.4826098960 


47.9016739726 


1.5809466820 


16 




52.6240518400 


51.0439078268 


1.5801529270 


17 




55.7655107528 


54. 1856175233 


1.5799006997 


18 




58.9069839241 


57.3278564233 


1.5791338198 


19 




62.0484691886 


60.4698329102 


1,5786416726 


20 




65.1899647992 


63.6115127150 


1.5784567241 
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TABLE III 

VALUES OF Z *** Jq(X) 



ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMU LAE 

p(jT(in = i+i/u( ( (a+3T(i n/2)? x(i + i) * a+tt/ \/p(xm) 



ERROR 


IS THE DIFFERENCE 


BETWEEN TABULATED 


AND COMPUTED ROOTS 


ROOT 

NO. 


TABULATED 

ROCT 


COMPUTED 

ROOT 


ERROR 


1 


5.5200781103 


5. 5217120022 


.0016338920 


2 


8.6537279126 


8.6538810786 


.0001531661 


3 


11.7915344390 


1 1. 7915695005 


.0000350617 


4 


14.9309177082 


14. 9309296827 


.0000119747 


5 


18.0710639674 


18. 0710691069 


.0000051396 


6 


21.2116366294 


21.2116391868 


.0000025579 


7 


24.3524715304 


24. 352472942B 


.0000014123 


8 


27.4934791317 


27.4934799736 


.0000008423 


9 


30.6346064676 


30. 6346070007 


.0000005335 


10 


33.7758202134 


33.7758205654 


. .0000003527 


1 1 


36.9170983527 


36.9170985967 


.0000002443 


12 


40.0584257636 


40.0584259368 


.0000001732 


13 


43.1997917127 


43. 1997918384 


.0000001260 


. 14 


46.3411883712 


46.3411884652 


.0000000946 


15 


49.4826098960 


49. 4826099686 


.0000000730 


16 


52.6240518400 


52.6240518950 


.0000000556 


17 


55.7655107528 


55. 7655107975 


.0000000452 


18 


58.9069839241 


58.9069839586 


.0000000349 


19 


62.0484691886 


62. 0484692166 


.0000000280 


20 


65.1899647992 


. 65.1899648216 


.0000000226 
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table IV 

VALUES OF $ *** Jq{X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTYt BASED ON THE RECURSION FORMULA 
p($(i)i = ui/4A)^n) ^(i + i) » u*r/ V?[ 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 



I- 



X( I ) ) 



ROOT 

NO* 


TABULATED 

ROOT 


COMPUTED 

ROOT 


ERROR 


1 


5.5200781103 


5.5172325156 


-.0028455946 


2 


8.6537279126 


8.6534819095 


-.0002460030 


3 


11.7915344390 


11.7914791510 


-.0000552879 


4 


14.9309177082 


14.9308989581 


-.0000187501 


5 


18.0710639674 


18. 0710559469 


-.0000080201 


6 


21.2116366294 


21. 2116326448 


-.0000039842 


7 


24.3524715304 


24.3524693325 


-.0000021976 


8 


27.49347913lt 


27.4934778223 


-.0000013094 


9 


30.6346064676 


30.6346056401 


-.0000008275 


10 


33.7758202134 


33.7758196630 


-.0000005496 


n 


36.9170983527 


36.9170979746 


-.0000003774 


12 


40.0584257636 


40.0584254945 


-.0000002689 


/ 13 


43.1997917127 


43.1997915152 


-.0000001969 


14 


46.341 1883712 


46.3411882240 


-.0000001468 


15 


49.4326098960 


49.4826097042 


-.0000001110 


16 


52.6240518400 


52.6240517525 


-.0000000871 


17 


55.7655107528 


55.7655106848 


-.0000000672 


18 


58.9069839241. 


58.9069838692 


-.0.000000548 


19 


62.0484691886 


62.0484691439 


-.0000000444 


20 


65.1899647992 


65. 1899647620 


-.0000000366 
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TABLE V 

VALUES OF in *•* Jq{X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

pixTm) * (P(x(i))+p($m))/2 Ttn+n- a+^ VpTxTm) 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 



ROOT 

NO. 


TABULATED 

ROOT 


COMPUTED 

ROOT 


ERROR 


1 


5.5200781103 


5.5254224626 


.0053443524 


2 


8.6537279126 


8.6549613329 


^0012334206 


3 


11.7915344390 


11.7919896499 


.0004552110 


u 


14.9309177082 


14.9311330232 


.0002153153 


5 


18.0710639674 


18.0711822445 


.0001182773 


6 


21.2116366294 


21.2117084200 


.0000717907 


7 


24.3524715304 


24.3525183247 


.0000467948 


8 


27.4934791317 


27.4935113089 


.0000321774 


9 


30.6346064676 


30. 6346295336 


.0000230662 


10 


33.7758202134 


33.7758373059 


.0000170928 


1 1 


36.9170983527 


36.9171113698 


.0000130179 


12 


40.0584257636 


40.0584359039 


.0000101405 


13 


43.1997917127 


43. 1997997640 


.0000080521 


14 


46.341 1883712 


46.3411948718 


.0000065007 


15 


49.4826098960 


49,4826152194 


.0000053242 


16 


52.6240518400 


52.6240562536 


.0000044135 


17 


55.765L‘107528 


55.7655144539 


.0000037015 


.18 


58.9069839241 


58.9069870561 


.0000031325 


19 


62.0484691886 


62.0484718634 


.0000026751 


20 


65.1899647992 


65. 1899671014 


.0000023024 
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TABLE VI 



VALUES OF X2 *** Jq(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 

ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

p(”x2(i)) =(p(7( I) )+P(x( I ) ) )/2 «(i + i) » A+'^/VpufTnT” 

ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 



ROOT 

NO. 


TABULATED 

ROOT 


COMPUTED 

ROOT 


ERROR 


1 


5.5200781103 


5.5194687883 


”.0006093219 


2 


8.6537279126 


8.6536814596 


-.0000464527 


3 


11.7915344390 


11.7915243239 


-.0000101150 


U 


14.9309177082 


14.9309143201 


-.0000033879 


5 


18.0710639674 


18. 0710625267 


-.0000014403 


6 


21.2116366294 


21.2116359160 


-.0000007132 


7 


24.3524715304 


24.3524711374 


-.0000003927 


8 


27.4934791317 


27. 4934788980 


-.0000002336 


9 


30.63*16064676 


30.6346063204 


-.0000001470 


10 


33.7758202134 


33.7758201 146 


-.0000000984 


1 1 


36.9170983527 


36.9170982856 


-.0000000665 


12 


40.0584257636 


40.0584257152 


-.0000000478 


13 


43.1997917127 


43. 1997916764 


-.0000000355 


14 


46.3411883712 


46. 3411883442 


-.0000000261 


15 


49.4826098960 


49.4826098764 


-.0000000190 


16 


52.6240518400 


52.6240518242 


-.0000000158 


17 


55.7655107528 


55.7655107416 


-.0000000110 


18 


58.9069839241 


58.91369839139 


-.0000000099 


19 


62.0484691886 


62.0484691802 


-.0000000082 


20 


65.1899647992 


65. 1899647918 


-.0000000070 
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TABLE VII 

CUMULATIVE ERRORS FOR 'x *** Jq(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER the first ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 



FROM 

FROM 


THE VALUE CF X 
THE FIRST ROOT. 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


ROOT 


TABULATED 


COMPUTED 


ERROR 


NO* 


ROCT 


ROOT 




1 


5.5200781103 


5. 5336719233 


.0135938131 


2 


8.6537279126 


8.6700534131 


*0163255008 


3 


11.7915344390 


1 1. 8088337616 


.0172993229 


4 


14.9309177082 


14.9486705526 


.0177528444 


5 


18.0710639674 


18. 0390637650 


«01 79997976 


6 


21.2116366294 


21.229785478? 


.0181488496 


7 


24.3524715304 


24. 3707171553 


.0182456251 


8 


27.4934791317 


27.51 1791 1 105 


.0183119791 


9 


30.6346064676 


30. 6529659065 


.0183594390 


10 


33.7753202134 


33.7942147609 


.0183945483 


11 


36.9170983527 


36. 9355196003 


.0184212479 


12 


40.0584257636 


40.0768677853 


.0184420225 


13 


43.1997917127 


43.2182502151 


.0184585024 


14 


46*3411063712 


46.3596601672 


.0184717962 


15 


49.4826098960 


49. 5010925708 


.0184826752 


16 


52.6240518400 


52.6425435282 


.0184916886 


17 


55.7655107528 


55. 7840099940 


.0184992421 


18 


58.9069839241 


58.9254895560 


.0185056320 


19 


62.0484691886 


62.0669802753 


.0185110874 


20 


65.1899647992 


65.2084805798 


.0185157806 
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TABLE VIII 

CUMULATIVE ERRORS FOR X *** Jq(X) 



ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 



AFTER 


THE FIRST ROOT 


IS COMPUTED, THE NEXT 


ROOT IS COMPUTED 


FROM 


THE VALUE OF X 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


FROM THE FIRST ROOT. 






ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


5.520078 1103 


5.5217120022 


.0016338920 


2 


8.6537279126 


8.6555185516 


.0017906391 


3 


1 1.791534L390 


1 1.7933614515 


.0018270126 


l4 


14.9309177082 


14.9327572957 


.0018395875 


5 


18.0710639674 


18. 0729090152 


.0018450481 


6 


21.2116366294 


21.2134844260 


.0018477967 


7 


24.3524715304 


24.3543208619 


.0018493315 


8 


27.4934791317 


27.4953293884 


.0018502570 


9 


30.6346064676 


30.6364573166 


.0018508493 


10 


33.7758202134 


33.7776714582 


.0018512453 


11 


36.9170983527 


36.9189498741 


.0018515220 


12 


40.0584257636 


40.0602774834 


.0018517202 


13 


43.1997917127 


43.2016435778 


.0018518660 


14 


46.341 1883712 


46.3430403462 


.0018519759 


15 


49.4826098960 


49.4844619567 


.0018520612 


16 


52.6240518400 


52. 6259039668 


.0018521272 


17 


55.7655107528 


55.7673629336 


.0018521811 


18 


58.9069839241 


58. 9088361468 


.0018522234 


19 


62.0484691886 


62.0503214458 


.0018522573 


20 


65.1899647992 


65.1918170843 , 


.0018522854 
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TABLE IX 

CUMULATIVE ERRORS FOR X *•* Jq(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 
FROM THE VALUE OF $ LAST OBTAINED. ERROR IS CUMULATIVE 
FROM THE FIRST ROOT. 



ROOT 

NO. 


TABULATED 

ROCT 


COMPUTED 

ROOT 


ERROR 


1 


5.5200781103 


5. 5172325156 


-.0028455946 


2 


8.6537279126 


8.6506294168 


-.0030984957 


3 


11.7915344390 


1 1.7883782729 


-.0031561660 


U 


14.9309177082 


14.9277417252 


-.0031759829 


5 


18.0710639674 


18.067879399S 


-.0031845679 


6 


21.2116366294 


21.2084477427 


-.0031888865 


7 


24.3524715304 


24.3492802321 


-.0031912980 


8 


27.4934791317 


27.4902863791 


-.0031927526 


9 


30.6346064676 


30. 6314127846 


-.0031936828 


10 


33.7758202134 


33.7726259045 


-.0031943080 


n 


36.9170983527 


36.9139036089 


-.0031947432 


12 


40.0584257636 


40.0552307060 


-.0031950568 


. 13 


43.1997917127 


43. 1965964232 


-.0031952894 


14 


46.3411883712 


46.3379929066 


-.0031954644 


15 


49.4826098960 


49, 4794142973 


-.0031955985 


- 


52.6240516400 


52.6208561352 


-.0031957048 


17 


55.7655107528 


55.7623149641 


-.0031957879 


16 


58.9069839241 


58.9037880665 


-.0031958566 


19 


62.0484691886 


62.0452732751 


-.0031959134 


20 


65.1899647992 


65. 1867688391 


-.0031959599 
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table X 

CUMULATIVE ERRORS FOR xT *** Jq(X) 



ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FOR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 



AFTER THE FIRST ROOT 


IS COMPUTED. THE NEXT 


ROOT IS COMPUTED 


FROM 


THE VALUE OF xi 


LAST obtained. ERROR 


IS CUMULATIVE 


FROM 


THE FIRST ROOT. 






ROOT 

NO. 


TABULATED 

ROOT 


COMPUTED 

ROOT 


ERROR 


1 


5.5200781103 


5.5254224626 


.0053443524 


.2 


8.6537279126 


8.6603153765 


.0065874640 


3 


11.7915344390 


11. 7985812158 


.0070467769 


4 


14.9309177082 


14.9381818194 


.0072641112 


5 


18.0710639674 


18.0784474835 


.0073835164 


6 • 


21.2116366294 


21. 2190926266 


.0074559972 


7 


24.3524715304 


24.3599747745 


.0075032441 


8 


27.4934791317 


27.5010148650 


.0075357334 


9 


30.6346064676 


30. 6421654909 


.0075590234 


10 


33.7758202134 


33.7833964955 


.0075762823 


1 1 


36.9170983527 


36. 9246877786 


.0075894268 


12 


40.0584257636 


40. 0660254285- 


.0075996651 


13 


43.1997917127 


43.2073995080 


.0076077954 


14 


46.3411883712 


46. 3488027304 


.0076143594 


15 


49.4826098960 


49.4902296308 


.0076197354 


16 


52.6240518400 


52.6316760313 


.0076241914 


17 


55.7655107528 


55. 7731386814 


.0076279288 


18 


58.9069839241 


58.9146150155 


.0076310916 


19 


62.0484691886 


62.0561029809 


.0076337926 


20 


65.1899647992 


65. 19760091 60 


.0076361171 
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TABLE XI 

CUMULATIVE ERRORS FOR 1(2 *** Jq ( X ) 



ZEROS OF THE BESSEL FUNCTION OF ORDER ZERO FQR ROOTS 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 



AFTER THE FIRST ROOT 


IS COMPUTED, THE NEXT 


ROOT IS COMPUTED 


FROM 


THE VALUE OF X2 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


FROM 


THE FIRST ROOT. 






ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


5.5200781103 


5.5194687883 


-.0006093219 


2 


8.6537279126 


8.6530707316 


-.0006571810 


3 


11.7915344390 


n. 7908666495 


-.0006677894 


4 


14.9309177082 


14.9302463082 


-.0006714000 


5 


18.0710639674 


18.0703910086 


-.0006729586 


6 


21.2116366294 


21. 2109628869 


-.0006737421 


7 


24.3524715304 


24.3517973502 


-.0006741801 


8 


27.4934791317 


27.4928046870 


-.0006744443 


9 


30.6346064676 


30.6339318538 


-.0006746134 


10 


33.7758202134 


33.7751454851 


-.0006747281 


1 1 


36.9170983527 


36.9164235452 


-.0006748069 


12 


40.0584257636 


40.0577508984 


-.0006748646 


13 


43.1997917127 


43. 1991168046 


-.0006749080 


14 


46.3411883712 


46. 3405134305 


-.0006749401 


15 


49.4826098960 


49.4819349311 


-.0006749644 


16 


■ 52.6240518400 


52.6233768547 


-.0006749846 


17 


55.7655107528 


55.7648357525 


-.0006749997 


10 


58.9069839241 


58.9063089108 


-.0006750130 


19 


62.0484691886 


62.0477941642 


-.0006750239 


20 


65.1899647992 


65. 1892897654 


-.0006750335 
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TABLE XII 

VALUES OF *X *** (X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 

ONE through twenty, BASED ON THE RECURSION FORMULAE 

p(x(D) = i+n-4N^)/u(x(i ))^ x(i + i) = a+ 71/ i/p(x( m 



ERROR 


IS THE DIFFERENCE 


BETWEEN TABULATED 


AND COMPUTED ROOTS 


ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


7.0155866698 


6.9976376109 


•“.0179490588 


2 


10.1734681348 


10. 1686351451 


-.0048329895 


3 


13.3236919360 


13.3217202593 


-.0019716767 


4 


16.4706300506 


16. 4696368496 


-.0009932009 


5 


19.6158585101 


19.6152890925 


-.0005694173 


6 


22.7600853802 


22. 7597279320 


-.0003574477 


7 


25.9036720875 


25.9034342724 


-.0002378150 


8 


29.0468205335 


29.0466620042 


-.0001665289 


9 


32.1896799095 


32. 1095587798 


-.0001211293 


10 


35.3323075492 


35. 3322166996 


T. 0000908489 


n 


38.4747662339 


38.4746963540 


-.0000698795 


12 


41.6170942122 


41. 6170393117 


-.0000549003 


13 


44.7593189972 


44.7592750816 


-.0000439151 


14 


47.9014608869 


47.9014252098 


-.0000356767 


15 


51.0435351832 


51.0435058065 


-.0000293767 


16 


54.1855536401 


54. 1855291631 


-.0000244770 


17 


57.3275254359 


57.3275048267 


-.0000206090 


10 


60.4694578443 


60. 4694403261 


-.0000175173 


19 


63.6113566970 


63.61 13416851 


-.0000150119 


20 


66.7532267310 


66.7532137688 


-.0000129620 
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table XIII 

VALUES OF 9^ *** Jj(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 



CHI IS THE VA LUE X MUST ASSUME TO CONVERGE TO B 
% - ( {B-A)/2) { >/ { l-4N^)/( (B-A ) 



ERROR 


IS THE DIFFERENCE 


BETWEEN X AND CHI 




ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROCT 


ROOT 




1 


7.0155866698 


5.3312873839 


1.6663502270 


2 


10.1734681348 


8. 5374782400 


1.6311569051 


3 


13.3236919360 


1 1.7071 112888 


1.6146089705 


4 


16.4706300506 


14. 8645464401 


1.6050904095 


5 


19.6158585101 


18.0163550223 


1.5989340702 


6 


22.7600853802 


21. 1610798943 


1.5986480376 


7 


25.9036720875 


24. 3119656080 


1.5914686644 


8 


29.0468285335 


27. 4576333771 


1.5890286271 


9 


32.1896799095 


30.6024408266 


1.5871179532 


10 


35.3323075492 


33. 7466444736 


1.5855722260 


n 


38.4747662339 


36. 8904200587 


1.5842762953 


12 


41.6170942122 


40. 0338375270 


1.5832017846 


13 


44.7593189972 


43. 1769980257 


1.5822770558 


14 


47.9014608869 


46. 3199358247 


1.5814893851 


15 


51.0435351832 


49.4627229534 


1.5807828531 


16 


54.1855536401 


52. 6053818800 


1.5801472832 


17 


57.3275254359 


55.7479370888 


1.5795677379 


18 


60.4694578443 


58. 890202 9432 


1.5792373829 


19 


63.6113566970 


62.0326487795 


1.5786929056 


20 


66.7532267310 


65. 1751772780 


1.5780364908 
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TABLE 

VALUES OF r 



ZEROS 


OF THE BESSEL FUNCTION 


ONE THROUGH TWENTY, BASED ON 
P(Xm) = 1 + (1-UN^)/4((A+X( I) 


ERROR 


IS THE DIFFERENCE 


BETWE 


ROOT 

NO. 


TABULATED 

ROOT 


r 

u 


1 


7.0155866698 


7. 


2 


10.1734681348 


10. 


3 


13.3236919360 


13. 


U 


16.4706300506 


16. 


5 


19.6158585101 


19. 


6 


22.7600853802 


22. 


7 


25.9036720875 


25. 


8 


29.0468285335 


29. 


9 


32.1896799095 


32. 


10 


35.3323075492 


35. 


1 1 


38.4747662339 


38. 


12 


41.6170942122 


41. 


13 


44.7593189972 


44. 


14 


'47.9014608869 


47. 


15 


51.0435351832 


51. 


16 


54.1855536401 


54. 


17 


57.3275254359 


57. 


18 


60.4694578443 


60. 


19 


63.6113566970 


63. 


20 


66.7532267310 


66. 



XIV 

*** Jj^(X) 

OF ORDER ONE FOR ROOTS NUMBER 

THE RECURSION FORMULAE 

)/2)^ 7ll + n = A+7Z/ /pixi I ) ) 

EN TABULATED AND COMPUTED ROOTS ' 



OMPUTED 

ROOT 


ERROR 


0141418219 


-.0014448478 


1732513588 


-.0002167758 


3236309085 


-.0000610273 


4706066190 


-.0000234314 


6158476649 


-.0000108451 


7600786900 


-.0000066901 


9036688204 


-.0000032669 


0468265270 


-.0000020061 


1896786084 


-.0000013006 


3323066691 


-.0000008800 


4747656174 


-.0000006158 


6170937680 


-.0000004442 


7593186684 


-.0000003282 


9014606383 


-.0000002479 


0435349923 


-.0000001903 


1855534911 


-.0000001483 


3275253186 


-.0000001170 


4694577483 


-.0000000957 


6113566197 


-.0000000766 


7532266676 


-.0000000615 
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table XV 

VALUES OF X *** J_^(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 

ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE^ 

p{5m> » 1 ♦ { 1 - 4N^)/ UAX( I) ^(I + l) » A +X/^P(^( I ) ) 
ERROR IS THE DIFFERENCE BETWEEN TABULATED AND COMPUTED ROOTS 



ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


7.0155866698 


7. 0180468041 


.0024601344 


2 


10.1734681348 


10. 1738162260 


.0003480912 


3 


13.3236919360 


13.3237883153 


.0000963795 


u 


16.4706300506 


16.4706667955 


.0000367450 


5 


19.6158535101 


19. 6158754556 


.0000169458 


6 


22.7600853802 


22.7600932522 


.0000078722 


7 


25.9036720875 


25.9036771739 


.0000050868 


8 


29.0468285335 


29.0468316572 


.0000031241 


9 


32.1896799095 


32. 1896819314 


.0000020222 


10 


35.3323075492 


35.3323089145 


.0000013660 


n 


38.4747662339 


38.4747671895 


.0000009564 


12 


41.6170942122 


41. 6170949005 


.0000006892 


13 


44.7593189972 


44.7593195057 


.0000005093 


14 


47.9014608869 


47.9014612706 


.0000003842 


15 


51.0433351832 


51.0435354784 


.0000002958 


16 


54.1855536401 


54. 1855538711 


.0000002316 


17 


57,3375254359 


57.3275256194 


.0000001642 


18 


60.4694576443 


60.4694579896 


.0000001462 


19 


• 63.6113566970 


63.6113568163 


.0000001 198 


20 


66.7532267310 


•66. 7532.268297 


,0000000997 



TABLE XVI 

VALUES OF ri *** J|(X) 



ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FqR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULA^ 

p(ri(i)) - ( p( x( I n+F (^( I ) ) )/ ? jTfd) * M+7C/V p(3(T( I) ) 



ERROR 


IS THE DIFFERENCE 


BETWEEN TABULATED 


AND COMPUTED ROOTS 


ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


7.0155866698 


7.0077901660 


-.0077965038 


2 


10.17346815411 


10. 1712216930 


-.0022464416 


3 


13.3236919360 


13. 3227535994 


-.0009383365 


4 


16.470630C506 


16.4701516442 


-.0004784064 


5 


19.6158585101 


19. 6155822147 


-.0002762953 


,6 


22.7600853802 


22.7599105686 


-.0001748114 


7 


25.9036720875 


25.9035557127 


-.0001163747 


8 


29.0468285335 


29.0467468258 


-.0000817076 


9 


32.1896799095 


32. 1896203524 


-.0000595564 


10 


35.3323075492 


35.3322628057 


-.0000447430 


1 1 


38.4747662339 


38.4747317713 


-.0000344625 


12 


41.6170942122 


41.6170671051 


-.0000271062 


13 


44.7593189972 


44.7592972936 


-.0000217033 


14 


47.9014608869 


47.9014432402 


-.0000176465 


15 


51.0435351832 


51.0435206424 


-.0000145406 


16 


54.1855536401 


54. 1855415171 


-.0000121228 


17 


57.3275254359 


57. 3275152231 


-.0000102125 


18 


60.4694578443 


60.4694491578 


-.0000086856 


19 


63.6113566970 


63. 6113492502 


-.0000074461 


20 


66.7532267310 


66.7532202993 


-.0000064312 



TABLE XVII 
VALUES OF ^2 ♦** 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

p(xTm) = (p(x( I >)+p(^(i n )/2 )^m = a+7T/ \/p c iTa n ) ) 



ERROR 


IS THE DIFFERENCE 


BETWEEN TABULATED 


AND COMPUTED ROOTS 


ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


7.0155866698 


7. 0160913781 


.0005047083 


2 


10.1734681348 


10. 1735337211 


.0000655864 


3 


13.3236919360 


13. 3237096060 


,0000176702 


4 


16.4706300506 


16. 4706367063 


.0000066559 


5 


19.6158585101 


19.6158615602 


.0000030501 


6 


22.7600853802 


22. 7600859711 


.0000005910 


7 


25.9036720875 


25. 9036729969 


,0000009099 


8 


29.0468285335 


29. 0468290923 


.0000005590 


9 


32.1896799095 


32. 1896802699 


.0000003608 


10 


35.3323075492 


35. 3323077913 


.0000002430 


1 1 


38.4747662339 


38.4747664034 


.0000001703 


12 


41.6170942122 


41. 6170943342 


.0000001225 


13 


44.7593189972 


44.7593190875 


.0000000905 


. 14 


47.9014608869 


47. 9014609549 


.0000000681 


15 


51.0435351832 


51.0435352353 


.0000000527 


16 


54.1855536401 


54. 1855536811 


.0000000416 


17 


57.3275254359 


57.3275254695 


.0000000336 


18 


60.4694578443 


60. 4694578694 


.0000000253 


19 


63.6113566970 


63.6113567185 


.0000000216 


20 


66.7532267310 


66. 7532267496 


.0000000191 
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TABLE XVIII 

CUMULATIVE ERRORS FOR 'J? »** J^tX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 



FROM 


THE VALUE OF 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


FROM 


THE FIRST RCOT. 






ROOT 

NO. 


TABULATED 

ROOT 


COMPUTED 

ROOT 


ERROR 


1 


7.0155866698 


6. 9976376109 


“.0179490588 


2 


10.1734681348 


10. 1507267661 


-.0227413687 


3 


13.3236919360 


13.2990017352 


-.0246902007 


4 


16.4706300506 


16. 4449597476 


-.0256703027 


5 


19.6158585101 


19.5896268408 


-.0262316693 


6 


22.7600853802 


22. 7335015247 


-.0265838551 


7 


25.9036720875 


25. 8768550311 


-.0268170560 


8 


29.0468285335 


29. 0198475332 


-.0269810002 


9 


32.1896799095 


32. 1625796901 


-.0271002193 


10 


35.3323075492 


35.3051179312 


-.0271896177 


11 


38.4747662339 


38.4475078629 


-.0272583706 


12 


41.6170942122 


41. 5897818329 


-.0273123789 


. 13 


44.7593189972 


44.7319634212 


-.0273555757 


14 


47.9014608869 


47. 8740702206 


-.0273906654 


15 


51.0455351832 


51. 0161 156254 


-.0274195571 


16 


54.1855536401 


54. 1581 1001 14 


-.0274436282 


17 


57.32?5254359 


57. 3000615416 


-.0274638942 


18 


60.4694578443 


60.4419767251 


-.0274811186 


19 


63.61 13566970 


63.5838608174 


-.0274958793 


20 


66.7532267310 


66.7257181071 


-.0275086237 
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table XIX 

CUMULATIVE ERRORS FOR X **• JiIX) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 



FROM 


THE VALUE OF X 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


FROM 


THE FIRST RCOT. 






ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 


1 


7.0155866698 


7.0141418219 


-.0014448478 


2 


10.1734681348 


10. 1718119476 


-.0016561872 


3 


13.3236919360 


13. 3219771462 


-.0017147896 


4 


16.4706300506 


16.4688930572 


-.0017369931 


5 


19.6158585101 


19.6141 113704 


-.0017471395 


6 


22.7600853802 


22.7583319838 


-.0017533960 


7 


25.9036720875 


25.90191671 11 


-.0017553763 


8 


29.0468285335 


29. 0450713504 


-.0017571828 


9 


32.1896799095 


32. 1879215697 


-.0017583393 


10 


35.3323075492 


35.3305484364 


-.0017591120 


1 1 


38.4747662339 


38. 4730065875 


-.0017596461 


12 


41.6170942122 


41.6153341858- 


-.0017600260 


13 


44.7593189972 


44.7575586932 


-.0017603031 


14 


47.9014608869 


47. 8997003762 


-.0017605102 


15 


51.0435351832 


51. 04177451 60 


-.0017606668 


16 


54. 1855536401 


54. 1837928528 


-.0017607870 


17 


57.32/5254359 


57. 3257645555 


-.0017608804 


IS 


60.4694578443 


60.4676968884 


-.0017609558 


19 


63.6113566970 


63. 6095956815 


-.0017610151 


20 


66.7532267310 


66.7514656689 


-.0017610620 



TABLE XX 

CUMULATIVE ERRORS FOR ^ *** Jj(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE through twenty, BASED ON THE RECURSION FORMULAE 



AFTER 


THE FIRST ROOT 


IS COMPUTED, THE NEXT 


ROOT IS COMPUTED 


FROM 


THE VALUE OF ^ 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


FROM THE FIRST ROOT. 






ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


7.0155866698 


7.0180468041 


.0024601344 


2 


10. 1734681348 


10. 1762664455 


.0027983109 


3 


13.3236919360 


13. 3265823792 


.0028904433 


U 


16.4706300506 


16. 4735551225 


.0029250719 


5 


19.6158585101 


19.6187993325 


.0029408227 


6 


22.7600853802 


22. 7630333365 


.0029479564 


7 


25.9036720875 


25. 9066256429 


.0029535557 


8 


29.0468285335 


29. 0497848745 


.0029563413 


9 


32. 1 896799095 


32. 1926380284 


.0029581190 


10 


35.3323075492 


35. 3352668518 


.0029593028 


1 1 


38.4747662339 


38.4777263533 


.0029601197 


12 


41.6170942122 


41. 6200549109 


.0029606996 


13 


44.7593189972 


44.7622801 177 


.0029611211 


14 


■'47.9014608869 


47.9044223214 


.0029614345 


15 


51.0435351832 


51.0464968551 


.0029616724 


16 


54.1855536401 


54. 1885154955 


.0029618555 


17 


57.3275254359 


57.3304874348 


.0029619992 


18 


60.4694578443 


60. 4724199548 


.0029621108 


19 


63.6113566970 


63.6143188979 


.0029622011 


20 


66.7532267310 


66. 7561890054 


.0029622754 
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TABLE XXI 

CUMULATIVE ERRORS FOR 7 \ *** J]^(X) 

ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FQR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 

AFTER THE FIRST ROOT IS COMPUTED, THE NEXT ROOT IS COMPUTED 



FROM 

FROM 


THE VALUE OF XI 
THE FIRST RCOT. 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 




1 


7.0155866698 


7.0077901660 


“.0077965038 


2 


10.1734681348 


10. 1634497307 


“.0100184040 


3 


13.3236919360 


13.3127478317 


-.0109441041 


U 


16.4706300506 


16.4592144503 


-.01 14155999 


5 


19.6158585101 


19.6041707369 


-.01 16877730 


6 


22.7600853802 


22.7482254347 


“.01 18599452 


7 


25.9036720875 


25.8916985528 


-.01 19735343 


8 


29.0468285335 


29.0347745535 


-.0120539797 


9 


32.1896799095 


32. 1775672976 


-.0121126117 


10 


35.3323075492 


35. 3201508904 


-.0121566580 


1 1 


38.4747662339 


38.4625756508 


-.0121905831 


12 


41.6170942122 


41. 6048769467 


-.0122172653 


13 


44.7593189972 


44.7470803680 


-.0122386289 


lU 


47.9014608869 


47. 8892048877 


-.0122559989 


15 


51.0435351832 


51.0312648714 


-.0122703116 


16 


54.1855536401 


54. 1732713953 


-.0122822442 


17 


57.3275254359 


57.3152331384 


-.0122922968 


18 


60.4694578443 


60.4571569972 


-.0123008466 


19 


63.6113566970 


63.5990485204 


-.0123081761 


20 


66.7532267310 


66.7409122232 


-.0123145066 
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TABLE XXII 

CUMULATIVE ERRORS FOR )(2 **• Jj(X> 



ZEROS OF THE BESSEL FUNCTION OF ORDER ONE FOR ROOTS NUMBER 
ONE THROUGH TWENTY, BASED ON THE RECURSION FORMULAE 



AFTER THE FIRST ROOT 


IS COMPUTED, THE NEXT 


ROOT IS COMPUTED 


FROM 


THE VALUE OF Zz 


LAST OBTAINED. ERROR 


IS CUMULATIVE 


FROM 


THE FIRST RCOT. 






ROOT 


TABULATED 


COMPUTED 


ERROR 


NO. 


ROOT 


ROOT 


1 


7.0155866698 


7. 0160913781 


.0005047083 


2 


10.1734681348 


10. 1740364628 


.0005683281 


3 


13.3236919360 


13. 3242770866 


.0005851508 


4 


16.4706300506 


16.4712214330 


.0005913828 


5 


19.6158585101 


1 9. 6 164 52702 L 


.0005941927 


6 


22.7600853802 


22. 7606000152 


.0005946353 


7 


25.9036720875 


25.9042685339 


.0005964467 


a 


29.0468285335 


29. 0474254708 


.0005969373 


9 


32.1596799095 


32. 1902771583 


.0005972489 


10 


35.3323075492 


35. 3329050038 


.0005974551 


1 1 


38.4747662339 


38.4753638301 


.0005975968 


12 


41.6170942122 


41. 6176919080 


.0005976967 


13 


44.7593189972 


44.7599167656 


.0005977688 


14 


47. 9 u. *.605 86 9 


47. 9020587085 


.0005978224 


15 


51. 0435351 832 


51.0441330457 


.0005978626 


16 


54.1855536401 


54. 1861515343 


.0005978945 


|7 


57.0275254359 


57.3281233553 


.0005979196 


18 


60.4694578443 


60.4700557813 


.0005979377 


19 


63.6115566970 


63.6119546490 


.0005979527 


20 


66.7532267310 


66.7538246959 


.0005979659 
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